Spin-Electric Coupling in Molecular Magnets 
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We study the triangular antiferromagnet Cua in external electric fields, using symmetry group 
arguments and a Hubbard model approach. We identify a spin-electric coupling caused by an 
interplay between spin exchange, spin-orbit interaction, and the chirahty of the underlying spin 
texture of the molecular magnet. This coupling allows for the electric control of the spin (qubit) 
states, e.g. by using an STM tip or a microwave cavity. We propose an experimental test for 
identifying molecular magnets exhibiting spin-electric effects. 
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Single-molecule mae nets (SMMs) [H have emerged as 
a fertile testing ground for investigating quantum effects 
at the nanoscale, such as tunneling of magnetization 
0, H, 3], or coherent charge transport 0, @,13|, or the 
decoherence and the transition from quantum to clas- 
sical behavior Q. SMMs with antiferromagnetic cou- 
pling between neighboring spins are especially promising 
for the encoding and manipulation of quantum informa- 



tion 
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|12| . for they act as effective two-level 
systems, while providing additional auxiliary states that 
can be exploited for performing quantum gates, even in 
the presence of untunable couplings between the qubits 
[1^. Intra- and inter- molecular couplings of SMMs can 
be engineered by molecular and supra-molecular chem- 
istry |14Lenabhng a bottom-up design of molecule-based 
devices |15| . 

While the properties of SMMs can be modfied dur- 
ing the synthesis, control on the time scales required 
for quantum information processing remains a challenge. 
The standard spin-control technique is electron spin 
resonance (ESR) driven by ac magnetic fields Bac{t) 
d, [12, [13, For manipulation on the time scale of 

1 ns (Rabi frequency ^ 10^ s~^) Bac should be of the 
order of 10~^ T, which, however, is difficult to achieve. 
The spatial resolution of 1 nm, required for addressing a 
single molecule, is also prohibitively small. At these spa- 
tial and temporal scales, the electric control is preferable, 
because strong electric fields can be applied to small re- 
gions by using, for example, STM tips 3, 13, 21 1. Also, 
the quantized electric field inside a microwave cavity can 
be used 22, 23, 24, [2^ to control single qubits and to 



induce coupling between them even if they are far apart. 

Here we identify and study an efficient spin-electric 
coupling mechanism in SMMs which is based on an in- 
terplay of spin exchange, spin-orbit interaction (SOI), 
and lack of inversion symmetry. Spin-electric effects in- 
duced solely by SOI [26| have been proposed 27| and ex- 
perimentally demonstrated [28| in quantum dots. How- 
ever, these SOI effects scale with the system size L as 
[irt - making them irrelevant for the much smaller SMMs. 
Thus, additional ingredients-such as broken symmetries- 



must be present in SMMs for an efficient coupling be- 
tween spin and applied electric field. 

In the following, we demonstrate the possibility of such 
spin-electric effects in SMMs by focusing on a specific ex- 
ample, namely an equilateral spin triangle, Cua [2^. In 
this SMM, the low energy states exhibit a chiral spin 
texture and, due to the absence of inversion symmetry, 
electric fields couple states of opposite chirahty. More- 
over, SOI couples the chirahty to the total spin, and thus 
an effective spin-electric interaction eventually emerges. 

Spin-electric coupling. The low-energy states of Cua 
can be described in terms of an effective spin Hamil- 
tonian. There, three spins-1/2 (one for each Cu^+ 
ion) are coupled through Heisenberg and Dzyaloshinski- 
Moriya interaction, accounting for (super-) exchange and 
spin orbit interaction, resp. [29|, 
3 3 
Hq Jii+iSi ■ Si+i -I- ^ Dii+i • Si X Si+i. (1) 

1=1 i=l 

The Daft symmetry of the triangle implies several rela- 
tions between the coupling constants [30|. Since Jpq ^ 



5K and ID^ 



0.5 K, the Heisenberg terms deter- 



mine the gross structure of the energy spectrum, and 
the Dzyaloshinski-Moriya terms the fine one. In particu- 
lar, due to the antiferromagnetic coupling {Jpq > 0), the 
ground state multiplct corresponds to total spin S ~ 1/2, 
and the gap to the first excited S — S/2 quadruplet is 
Ah = 3 J/2. The S* = 1/2 subspace is spanned by the 
symmetry-adapted states \x,M), i.e., 

|±l,-fl/2) = (UTT>+e±ITiT>+eTlTti»/V3, (2) 

I ±1,-1/2) = {\m) + e±\m) + e^\im/V3,{3) 

with e± — Q^^^'^/^^ that are simultaneous eigenstates of 
the chirahty operator Cz and of Sz (total spin), for the re- 
spective eigenvalues x s-nd M. Here, we have introduced 
the chirahty vector C with components 

= (-2/3) (si • S2 - 2s2 • sa + sa • si), (4) 

Cj, = (2/V3)(si -82-83 -si), (5) 

Cz = (4/V3)8i • (82 X 8a). (6) 
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FIG. 1: (Color online) (a) Cus-triangle exposed to an electric 
field E created by, e.g., an STM-tip. For E = 0, the exchange 
couplings in the molecule are equal (light triangle). A finite 
E affects the (super-) exchange coupling in a directional way 
(dark triangle), (b, c) Low-energy S — 1/2 states of Cus in a 
magnetic field B, with the zero-field SOI splitting Aso=l K- 
Light (red) and dark (blue) lines correspond to states with 
X =-|-l( — 1). If B II 2 (b), the transitions induced by E (thin 
arrows) conserve for B jj" 2: (c), these transitions result in 
a change of spin orientation (thick arrows). 



They satisfy [Cfc,C/] = i2efe;„C„ and [Cfe, Si] = 0, and act 
as Pauli matrices in the \x = ±1) bases. The states 
Ix, M) are spht by the SOI (see below) into two Kramers 
doublets, I ± l,±l/2) and | ± l,Tl/2). 

Next, we study the effect of an electric field E on 
the Cu3-spins using general symmetry group arguments. 
The low-energy molecular states form '^E' and ir- 
reducible representations (IRs) of D^h, with the S = 
1/2-states lowest in energy and transforming as two 
E' IRs, while the higher energy S — 3/2-states span 
the four A'^ IRs [H, Is^]. The chiral states, Eqs. (pt5|l . 
form the standard bases of the i5'-representations, with 
\E'^,S^) = Ix ±1,'S'.,). An electric field E couples 
to Cu3 via eE • R, where e is the electron charge, and 
R = transforms according to the reducible vec- 

tor representation Dv = A'2 ® E' . The A'2 IR acts on 
Z and the two-dimensional IR E' acts on the (X,Y)- 
components in the Cu3-plane. The only non-zero matrix 



elements of R in our basis are e{E'_^^ Sz\X-\E'_, Sz) = 
e{E'_,Sz\X+\E'^,Sz) = 2id, with d real denoting the 
electric dipole coupling, and X± — ±X + lY being the 
irreducible tensor components of R. The zero matrix el- 
ements imply that the S — 3/2-states are unaffected in 
first order in E, and E^ has no effect on the spins at all. 
The resulting coupling between the E-field and chirality 
is described by 5H-E, = dE' ■ C\\, where E' = 7?.z((/»)E is 
rotated by ^ = 77r/6 — 29 about z, and C\\ = {Cx,Cy,Q). 

To emphasize that the spin-electric effect derived above 
is based on exchange, we reinterpret our results in terms 
of spin interactions. In an equilateral triangle, and in the 
absence of electric field, the spin Hamiltonian is given by 
Eq. ([1]) with equal exchange couplings Ji,,;+i = J. Using 
then Eqs. @ and (0), we find 

SHj, = ^ ^sin[2(l - z)^/3 + e]s, ■ s,+i, (7) 

^=l 

where 6 is the angle between an in-plane E-field and the 
vector ri2 pointing from site 1 to 2. This form of 5H^ 
shows that the E-field lowers the symmetry by introduc- 
ing direction-dependent corrections to the exchange cou- 
plings Jii+i. E.g., if 6* = 7r/2, 5J2Z = <5J3i ^ (5Ji2. We 
mention that the eigenstates can be labeled by a partial 
spin sum quantum number S12 = 0, 1 of (si + 82)^. 

While the specific form of d in Eq. ([T]) depends on mi- 
croscopic details, the structure of SH-^ relies solely on 
symmetry arguments and the lack of inversion symme- 
try is crucial for the spin-electric coupling. Indeed, in 
inversion-symmetric SMMs, the electric field, which is 
odd under inversion, and the spin, which is even, can 
only couple through transitions between orbitals of op- 
posite parity. When the orbital state is quenched, as in 
the spin Hamiltonian description, there are no orbitals of 
opposite parity. Thus, in the symmetric case, there is no 
spin-electric effect up to linear order in E. 

Next, we turn to the SOI. The most general 
form of SOI allowed by the Dj,h symmetry reads, 
^so = ELi[4o^a^'4 + A^o(^^B^S- + Te'1s%)1 where 
Ago('^so) effective SOI coupling constant for the 

A2- (E'J.-) irreducible representation, and T^'^ (Te'i) is 
the corresponding irreducible tensor operator [s^l- Us- 
ing again symmetry group arguments, we find that the 
SOI Hamiltonian acting in the S" = 1/2 subspace reads 

5Hso = ^soCzSz, where Ago = a|o. Using Eq. ([6]), 
this general form can be reduced to the Dzyaloshinski- 
Moriya interaction given in Eq. ([1]) . 

The coupling to a magnetic field B is given by B • g • S, 
with the Bohr magneton absorbed in the g-tensor. Due to 
the Ds/j-symmetry, g is diagonal with components gxx = 
9yy ~ 5-L the Cus-plane and g^z = g\\ normal to it. 

Combining 5H^ and SHgo, we finally obtain the effec- 
tive low-energy Hamiltonian in the presence of SOI and 
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FIG. 2: (Color online) Low-energy spectrum of the Cua 
molecule, (a) Energy levels of the Cua molecule in an in- 
plane magnetic field By (black solid line), split as a static 
in-plane electric field is turned on (dashed red line and 
dotted blue line), (b) The electric dipole coupling d can be 
determined from the slope of ESR lines in a constant magnetic 
field [vertical lines in (a)] as a function of . 



electric and magnetic fields, 
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-dE'-C||. (8) 



From this we see that an in-plane E-field causes rota- 
tions of the chirahty pseudospin. To illustrate the role 
of B, we focus on the case E || rai, corresponding to 
5H-E, = —dECx ■ For B || z, the eigenstates coincide with 
those of Sz, and thus E will not induce transitions be- 
tween |±, 1/2) and |±, — 1/2), but will do so in subspaces 
of given M, see Fig.[TJ For B || z, instead, the system 
eigenstates for E = are no longer eigenstates of Sz, and 
thus the electric field induced transitions result in spin 
flips, see Fig.[T]Jc). 

We emphasize that the electric dipole coupling d in 
Eq. ^ can be directly accessed in experiments, e.g., 
by standard ESR measurements [ll] in static electric 
fields, see Fig. [5] While a microscopic evaluation of 
d requires an ab-initio approach, which is beyond the 
scope of this work, we can estimate d, |E| and the spin- 
manipulation (Rabi) time resulting from Eq. ([5]) as fol- 
lows. For d between c?,„in = \Q~^eRi2 and dmax — eRi2 
and for E ^ lO^kV/cm, obtainable near an STM tip, 
the Rabi time is TRabi ~ 0.1 — 10'^ ps. The condition 
dE <C Ah for the validity of -ff^g " in Eq. ^ provides an- 
other lower bound on the spin-manipulation time, namely 

mill ^ in™ 
^Rabi ~ J^LI pi'- 

Hubbard approach. In order to gain further in- 
sight into the interplay between the exchange 
interaction and the electric field, we introduce 
an iVs-site Hubbard model of the triangular spin 
chain. The corresponding Hamiltonian reads, Hn = 



h. c. 



where Ui is the repulsion on site ?, ta^i the hopping 



FIG. 3: (Color online) (a) Hubbard models A and B of the 
spin triangle. Model A: Ne = Ns=3, ti=t, Ui = U, and ei = 
al'i+i=0. Model B: N^^9, Ns^6, al,+^=0, a^+i^a, 
ti^t, e3fc-2 -esfc'-i = U3k-2-U3k'-i = U {k,k' = 1,2,3). 

(b) Overlap between the projected ground state of Hh + He 
(I'^i)), and the eigenstates of si • S2 corresponding to Si2=0 
(squares) and Si2 = l (triangles), as function of the angle 9 
between the triangle side 1-2 and an in-plane E. The filled 
(empty) symbols correspond to the A (B) model, whereas the 
dotted lines give the components of the SHe ground state. 
In both models, t/U = 0.1, eRE/U = 2.5 x 10"^, and a = 0.1. 

(c) Dependence of A21 on the amplitude E, for E || y and 
eREmax/U=2.5 x 10"^. Filled (empty) symbols refer to the 
A (B) model, squares (triangles) correspond to a=0.1 (q=0) 
two-site contributions. 



matrix element, a — tiii ^J^d X)i a- "-f ~ The 
coupling of the system to the electric field E is 

i/E = eE • ^ [n,,<,r, + (?»;4-i4<,c,+i,<, + h. c.)] . (9) 

In the single-site terms, the expectation value of the elec- 
tron position r in the Wannier state \(t>i) is identified 
with the corresponding ion position, {(j)i\r\(j)i) ~ r^. The 
two-site terms describe the electric-field assisted hop- 
ping of electrons between neighboring sites, with ru^i = 
{(f)i\r\(j>i+i) = a^l^_^_^rii+l+a-'-ezXr^i+l. We now focus on 
the two main mechanisms giving antiferromagnetic cou- 
pling, namely direct exchange and superexchange (mod- 
els A and B, Fig. 2(a)). In both cases, the low-energy 
subspace {Sq) is defined by the states (|q;°)) where the 
magnetic ions at the triangle vertices are singly occu- 
pied. For E = 0, the projection of these states onto Sq 
(l*i-8>) coincides with the 5" = 1/2 and S = 3/2 eigen- 
states of the Heisenberg Hamiltonian. The degeneracy in 
the iS* = 1/2 multiplet is lifted by an electric field. 

In Fig.[2{b), we show the overlap between the projected 
ground state 1*5) and the \S = 1/2, S'12 = 0, 1) states for 
a given Sz as function of the direction of E (angle 9). 
The results coincide with the ones from Eq. ([7]) , for both 
models A and B. In addition, we find that the splitting 
(A21 = E2 — El) between the two lowest energies varies 
by less than 5% with 9, in agreement with H^"^ that 
predicts no ^-dependence at all. 

In Fig.[3Ic) we isolate the contribution to A21 arising 
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from the single- and two-site terms. All these contribu- 
tions scale linearly with E for every 0. The dependence 
of A21 on t, however, is model-dependent. In particular, 
in model A, the contributions to A21 arising from the 
single- and two-site terms scale like {t/U)^ and {t/U), 
resp. Analogous power-law dependences are found in 
model B, where the single-site (two-site) contribution 
scales as {t/Uy {{t/U)^), and two-site terms dominate in 
both models when t <^U . Additional mechanisms, such 
as the relative displacements of the ions, can contribute 
to the coupling between spin and electric field. 

Spin coupling to cavity electric fields. Exchange cou- 
pling of SMMs has been demonstrated in dimers 31 1. 
The use of this short-range and (so far) untunable inter- 
action requires additional resources for quantum informa- 
tion processing Efficient spin-electric interaction, 
on the other hand, provides a route to long-range and 
switchable coupling between SMM qubits. In particular, 
microwave cavities are suitable for reaching the strong- 
coupling regime for various qubit systems 22, 2^, 24, 25j . 
Here, we propose to use such cavities to control single 
SMMs and, moreover, to couple the spin-qubits of dis- 
tant SMMs placed inside the same cavity. 

The interaction of a single SMM with the cavity field 
reads, SH-^ — c?Eq • C||(fe'l' + 6), where Eg is the ro- 
tated electric field of amplitude |Eo| oc ^Jhio/V inside 
the cavity of volume V 23], and b is the annihilation 
operator for the photon mode of frequency lo. In the 
rotating wave approximation, the low-energy Hamilto- 
nian of N SMMs interacting with the cavity mode is 
i/,_ph = H^^^ + ^ where 

H^^^ = AsoCp^S'p) H- B • ^ S(^) + dEoie'^P^ b^C^l'^ + h. c), 

(10) 

y and ipj = 9j + 77r/6. Here, J?s-ph 
reduces to the well-known Tavis-Cummings model j32| 
when the spins are in eigenstates of sf', and B || z. 
However, if B 2; it is possible to couple both chirali- 
ties and total spins of distant molecules. Typically, the 
electric fields in cavities are weaker, |Eo| ~ IV/cm for 
hw « 0.1 meV 
giving TRabi ~ 
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than the ones near STM tips, thus 
0.01 — 100 /iS. Obviously, decreasing the 
cavity volume V would give shorter TRabi- Coupling of 
distant SMMs can be controlled by tuning two given 
molecules in and out of resonance with the cavity mode, 
e.g., by applying additional local electric fields. Further 
effects such as the state transfer between stationary and 
flying qubits, or the SMM-photon entanglement, can be 
observed in a system described by Hg-ph- 

In conclusion, we find an exchange-based mechanism 
that couples electric fields to spins in triangular molecular 
antiferromagnets. While our results are derived for Cua, 
analogous symmetry arguments are expected to apply to 
other molecular ma gnet s that lack inversion symmetry, 
such as ¥15(13, Coalil, Dygdl, Mni2@,i,[3, etc. 
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